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ABSTRACT 

We present and discuss the results of new multiple-scattering radiative transfer calculations for three 
representative types of galactic environments, filled with either homogeneous or two-phase clumpy dust 
distributions. Extinction and scattering properties for two types of interstellar dust, similar to those 
found in the average diffuse medium of the Milky Way Galaxy (MW) and the Bar of the Small Mag- 
ellanic Cloud (SMC), are considered. The wavelength coverage extends from 1000 A to 30,000 A, with 
particular emphasis on the rest-frame UV. This makes these models especially applicable to starburst 
galaxies and Lyman-break galaxy samples. The examination of the models concentrates on the study 
of UV/visual/near-IR reddening effects, the wavelength dependence of attenuation, and on the changes 
that arise from the transition from homogeneous to clumpy dust distributions in different star/dust 
geometries. Embedded dust, especially when clumpy, leads to saturation at fairly low reddening values 
with correspondingly gray attenuation functions. This makes the assessment of the attenuation of the 
far-UV flux from starburst galaxies difhcult, if only UV/visual/near-IR data are available. Existing data 
for UV-selected starburst galaxies indicate a range of UV attenuation factors of 0-150. Our models 
reproduce the "Calzetti Attenuation Law" , provided one adopts SMC-type dust and a clumpy shell-type 
dust distribution surrounding the starbursts. The average far-UV attenuation factor for the Calzetti 
sample is found to be 7.4. The only relatively reliable measure for the UV attenuation factor for single 
galaxies was found in the ratio of the integrated far-IR flux to the far-UV flux, measured near 1600 A, 
requiring the measurement of the entire spectral energy distribution of galaxies. 

Subject headings: dust extinction - radiation transfer 



1. INTRODUCTION 



Thronso n, & Capuano 1 992, h ereafter WTCj ; Bruzual 



The presence of dust in galaxies has a profound im 



al. 1988). According to WTC, a galaxy may be looked 



pact on observable characteristics such as their apparent 
UV/optical/near-IR luminosity, their spectral energy dis- 
tribution, their recognizable level of star-formation activ- 
ity, and their morphology. In addition, dust directly affects 
and often controls the physical conditions within the inter- 
stellar medium of galaxies and plays a major role in inter- 
stellar chemistry by providing the sites for the formation 
of interstellar molecular hydrogen and by regulating the 
interstellar radiation field. Quantitative estimates of the 
dust mass responsible for the observed effects are needed 
in order to yield gas-to-dust ratios in galaxies of differ- 
ent types, essential if we are to understand the chemical 
evolution of galaxies over cosmological time scales. 

Radiative transfer simulations often provide the most 
direct way to arrive at a quantitative assessment of the 
ways that lead to observed dust effects, provided such sim- 
ulations are sufficiently realistic to apply to actual physi- 
cal systems. All such simulations involve approximations. 
Hiatorically, thcac have ra nged from the a doption of ab 



upon as consisting of environments which are characterized 
by certain geometrical relations between the attenuating 
dust and the illuminating stars. For example, starburst 
regions may be approximated by regions occupied by lu- 
minous stars, surrounded on the outside by the remnants 
of molecular clouds from which the stars originated. Qui- 
escent parts of dusty disk galaxies may resemble systems 
in which stars and dust are fairly uniformly mixed, or, if 
one considers the same case at longer wavelengths where 
the dominant sources might be K-giants, a case where the 
scale height of sources is much larger than the scale height 
of the dust which is concentrated toward the center or cen- 



tral plane. As shown by WTC , equal amounts of dust in 



these different configurations produce very different red- 
dening and attenuation effects. For the analysis of effects 
resulting from changes in the inclination of disk galaxies, 
disk models with embedded dust have been developed, ei- 
ther in plane-parallel or doubly exponential geometry (e.g. 



rbing screen goomotrios (Holmbcrg 1958, 1975) in anal 



[Bruzual et al. 1988 ; Byun et al. 1994 ; di Bartolomeo et al ^ 
1996; Corradi et.al. 1996| ). All of the models mentioned so 



ogy to the study of extinction of starlight in our Galaxy, 
to the consideration of various geometries with embedded 
dust, including the effects of multiple scattering (Witt, 



far consistently illustrate the fact that a certain quantity of 
embedded dust produces much smaller attenuation effects 
than a foreground screen made from the same amount of 
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dust, and that the inclusion of multiple scattering reduces 
these effects even further by returning the scattered radi- 
ation to the optical radiation field, albeit with an altered 
spect ral and directional distribution. 



An approximation in common among all the models 1997 



presented in the present paper are particularly well suited 
to analyze the reddening effects and the corresponding 
UV attenuation in Lyman- brcak galaxies and galaxies seen 
in th e Hubble deep fields ([Dickinson 1997 ; [Mcurer ct al 



sumed dust distributions, which are either constant or of 
a monotonically varying continuous variety. In the first pa- 
per in the present series (Witt & Gordon 1996, hereafter 
WG96| ), we conducted an extensive exploration of the ef- 



^ 

mentimcd above is the homogeneous nature of their as- 1997; [Sawicki fc Yee 1998| , [Steidcl et al. 1999 ). They are 



■Icurer, Heckman, fc Calzctti 19991; jawicki ct al 



fects of dumpiness on the transport of radiation through 
a spherical system containing scattering dust and a single 
central source. This study, together with the res ults of 
earlie r work on transfer throug h clumpy media (e.g. Natta 
& Pa jiagia 1981 ; |Boissc 199oi ; [Hobson fc Scheucr 1993^ 



clearly revealed that the clumpy multi-phase structure of 
the interstellar medium (ISM) leads to a significant re- 
duction in the opacity produced by a giv en distr ibution 
and mass of interstellar dust. In addition, WG96 showed 



that the relationship between the effective optical depth of 
a clumpy distribution and that of an equivalent homoge- 
neous distribution is non-linear; this leads to an effective 
attenuation as a function of wavelength that differs from 
the opacity law intrinsic to the dust in the system. These 
conclusions have been further amplified in a recent investi- 
gation of the radiative transfer in the clumpy environment 
of young stellar objects by Wolf, Fischer, & Pfau (1998) 
and general investigation of clumpy radiative transfer by 
Varosi fc Dwek (1999). Since dumpiness is a reality of 
the ISM wherever it can be observed with sufficient spa- 
tial resolution, we are introducing it into the transfer of 
radiation in such galactic environments as investigated in 



WTC now. A preliminary version of some aspects of this 
work appeared in Gordon, Calzetti, & Witt (1997). 

In addition to adding models with a clumpy two-phase 
structur e, we have improved several other aspects of the 
original WTC study. Our knowledge of dust scattering 



properties in the UV has greatly improved since 1992, and 
with some confidence we are now modeling the transfer for 
15 (out of a total of 25) bandpasses in the wavelength range 



from 1000 A to 3000 A, compared to two in |WTC| . This 
improved UV coverage makes the current models particu- 
larly well suited for the analysis of starburst galaxies, both 
nearby and at high redshift, since these systems emit much 
of their stellar energy at UV wavelengths. Furthermore, 
most previous mult i- wavelength radiative transfer studies 
have assumed opacity functions and scattering properties 
representative of average Milky Way (MW) dust. Even 
within the MW dust c haracteristics var y from one environ- 



ment to another (e.g. Fitzpatrick 1999) and already in the 



opacity laws (Bianchi et al. 1996; 


Gordon fc Clayton 1998 


Misselt, Clayton, & Gordon 1999 


). In starburst galaxies 


in particular, it has been shown (Gordon ct al. 1997) that 



the opacity function intrinsic to the dust present there is 
much more similar to dust found typically in the Bar of 
the Small Magellanic Cloud (SMC). We therefore model 
all our galactic environments with both the average MW 
and SMC bar dust. Other radiative transfer studies which 
have included both MW and SMC bar type dust are those 
detailed by Gordon et al. (1997), Takagi, Arimoto, fc Van- 
sevicius (1999), and Ferrara et al. (1999). The models 



equally well usable for studying opacity effects in nearby 
dusty galaxies (e.g. Bcrlind et al. 1997; Gordon et al. 1997 



Kuchinskiet al. 1998| ; [Smith, Herter, fc Haynes 1998] ). The 
principal results from the current work are an extensive set 
of tables, accessible in the on-line version of this issue of 
The Astrophysical Journal or available by contacting the 
authors. We will present the details of the models in the 
following Section 2; selected results will be summarized in 
graphical form in Section 3; we will discuss the implica- 
tions of these results for currently ongoing galaxy studies 
in Section 4, to be followed by a summary in Section 5. 

2. MODEL 

In order to model the dust radiative transfer in com- 
plex systems such as galaxies, we use a radiative transfer 
model (nam ed DIRTY) based on Monte Carlo techniques 
( [Witt 1977 ). Use of such a model allows for arbitrary 
distributions of both the emitting sources (stars) and ab- 
sorbing and scattering material (dust). The full descrip- 
tion of DIRTY can be found in Gordon et al. (1999b). 
The inputs to DIRTY are the stellar distribution, the dust 
distribution, and the dust grain properties. With Monte 
Carlo techniques, individual photons are emitted isotrop- 
ically from stars according to the stellar distribution and 
followed through a dust distribution. Each photon's inter- 
action with the dust is parameterized by the dust optical 
depth (r), albedo, and scattering phase function. The 
optical depth determines where the photon interacts, the 
albedo gives the probability the photon is scattered from 
a dust grain, and the scattering phase function gives the 
angle at which the photon scatters. 

2.1. Dust and Stellar Distributions 



As was done in WTC, we have chosen to model spherical 
galactic environments instead of entire galaxies. A galac- 
tic environment describes the distribution of the stars and 
gas. In this paper we concentrate on three galactic envi- 
ronments: CLOUDY, DUSTY, and SHELL. These three 
geometries span the range of star/dust distributions and 
are pictorially represented in Figure The CLOUDY ge- 
ometry has dust extending to 0.69 of the system radius 
and stars extending to the system radius. The DUSTY 
geometry has both dust and stars extending to the system 
radius. This geometry represents a uniform mixture of 
stars and dust. The SHELL geometry has stars extending 
to 0.3 of the system radius and dust extending from 0.3 to 
1 of the system radius. This geometry was named nuclei 



in WTC 



These three geometries characterize the global star and 
dust distribution. The local distribution of the dust can be 
either homogeneous or clumpy. For this paper, w e define 
the clumpy distribution (using the definitions of WG96) 
to have a filling factor (jf) of 0.15, a low-density to high- 
density ratio (fe/fci) of 0.01, and a system size character- 
ized by a cube divided into 30'^ cubical bins {N = 30). 
A forthcoming third paper in this series will present the 
radiative transfer through doubly exponential disk/bulge 
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galaxies with clumpy ISM. 



Fig. 1. — Pictorial representation of the three geometries used 
in this paper (SHELL, DUSTY, and CLOUDY). The star symbols 
trace the stellar distribution and the grey solid regions give the dust 
distribution. 

2.2. Dust Grain Characteristics 

We have chosen to compute models using both Milky 
Way (MW) and Small Magellanic Cloud (SMC) dust. By 
using these two different types of dust, we are spanning the 
range of known dust extinction curves. For Monte Carlo 
radiative transfer, we need to know the dependence of the 
dust's optical depth (r), albedo (a), and scatterin g phase 
function on waveleng th. In this work, we used the Henyey 
& Greenstein (1941) phase function which is parameter- 
ized by the phase function weighted average of the cosine 
of the scattering angle {g — {cosO)). 

For the MW, we used dust which is associated with the 
diffuse interstellar medium and has Ry = Ay / E{B — V) = 
3.1. The wavelength dependence of r was determined from 
the parameterization of Milky Way extinctio n curv es by 
CardeUi, Clayton, & Mathis (1991, hereafter |CC^. The 
a and g values were determined by comparing empirical 
determinations of these numbers to a dust grain model 
which reproduces the R y = 3.1 extinction curve ( |KimJ 
Martin, & Hendry 1994). Figure ||a shows the empirical 
determinations, dust grain model values, and the numbers 
adopted in this paper. In determining the albedo values 
to use for this paper, we visually averaged the empirical 
determinations and then shifted the average to agree with 
the dust grain model values. The shift was done because 



most of the empirical determinations were done for reflec- 
tion nebulae which usually have higher values of Ry than 
3.1. The dust grain model predicts decreasing albedo's for 
decreasing i?y's. This is also consistent with recent work 




on the ultraviol et diffuse galactic light (Witt, Friedmann 



& Sasseen 1997). When we averaged the albedo values, we 
also made sure that Qsca {— o-Qoxt) was smooth across the 



2175 A bump as it is an absorption feature ( Calzetti et al 
plots Qoxt 



Figure 



Qsca, and Qabs for the MW 




albedo values used in this paper. The g values shown in 
Fig. 2b were taken as averages of the empirical determina- 
tions, as the g values from the ultraviolet diffuse galactic 
light work are not signiflcantly different from those of work 
on reflection nebulae. Also, the dust grain model still has 
problems in that it predicts the 2175 A bump as partly 
a scattering feature, contrary to observational evidence. 
The r, a, and g values adopted for MW dust in this paper 
are given m Table |l|. 

For the SMC, we used dust w hich is associated with th e 
star forming Bar of the SMC ([Gordon fc Clayton 1998| ). 
This dust produces extinction curves which arc quite dif- 
ferent from dust extinction c urves in the MW (CCM) or 
the Large Magellanic Cloud (Misselt, Clayton, fc Gordon 



1999]) . It lacks a 2175 A bump and is roughly linear with 
A^^. This type of dus t seems to be associated with intense 



star forming regions ( Gordon et al. 1997 Gordon fc Clay- 



iviissclt, Clayton, fc Gordon 1999|) . The r values 
were taken from an aver age of the three known ex tinction 
curves in the SMC Bar ( [Gordon fc Clayton 19*9^ ). There 
is no empirical work on the a and g values in the SMC. As 
there are also still problems with dust grain models repro- 
ducing the MW a and g values, we have chosen to start 
with the MW a values and then modify them to reflect the 
lack of a 2175 A bump. This time, we required that the 
Qsca and Qabs values did not show any features across the 
2175 A region. Figure ||b displays Qcxt, Qsca, and Qabs for 
the SMC albedo values used in this paper. The t, a, and 
g values adopted for SMC dust in this paper are given in 
Table |l|. 

3. RESULTS 

3.1. Structure of the Data Tables 

The results of our radiative transfer calculations are con- 
tained in Tables 2, 3, and 4, accessible in full in the online 
edition of this issue or available directly from the authors. 
Table 2 contains the data for the DUSTY galactic environ- 
ment, Table 3 those for the SHELL galactic environment, 
and Table 4 those for the CLOUDY environment. Only 
the beginning of Table 2 is shown in this paper in order 
to illustrate the type of data available in the full set of 
tables. Each model is identifled by its geometry (DUSTY, 
SHELL, CLOUDY), the intrinsic opacity function (MW, 
SMC), the radial extinction optical depth from the center 
to the edge of the dust environment at V {ry = 0.25, 0.5, 
...,50), assuming a constant density homogeneous distribu- 
tion, and the type of structure (homogeneous, clumpy) in 
the dust region. Each individual model contains calcula- 
tions for 25 wavelengths, ranging from 1000 A to 30000 A, 
shown in column 1. Column 2 identifles the optical depths 
for extinction at each wavelength, derived from the appli- 
cable intrinsic opacity function. Columns 3, 4, 5, and 6 
list, respectively, the attenuation optical depth (Tatt), the 
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Fig. 2. — The albedo (a) and scattering phase function asymmetry (g) values are plotted as a function of wavelength in (a) and (b), 
respectively. This plot includes values from empirical determinations from the literature, a dust grain model (Kim et al. 1994), and numbers 
adopted for this paper. 
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Fig. 3. — The various Q values for the MW and SMC are plotted in (a) and (b), respectively. 
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Table 1 
Dust Physical Properties 



MW SMC 





t/tv 




CL 




CP 

6 


t/tv 




d 




cr 
6 


1000.0 


5, 


.238 


0, 


.320 


0, 


,800 


9, 


,675 


0, 


,400 


0, 


,800 


1142.9 


3, 


.918 


0, 


,409 


0, 


,783 


7, 


,440 


0, 


,449 


0, 


,783 


1285.7 


3 


.182 


0, 


,481 


0, 


,767 


6, 


,068 


0, 


,473 


0, 


,767 


1428.6 


2. 


.780 


0, 


,526 


0, 


,756 


5, 


,167 


0, 


,494 


0, 


,756 


1571.4 


2, 


.584 


0, 


,542 


0, 


,745 


4, 


,536 





,508 


0, 


,745 


1714.3 


2, 


.509 


0, 


,536 


0, 


,736 


4, 


,074 


0, 


,524 


0, 


,736 


1857.1 


2, 


.561 


0, 


,503 


0, 


,727 


3, 


,700 


0, 


,529 


0, 


,727 


2000.0 


2, 


.843 


0, 


,432 


0, 


,720 


3 


,379 


0, 


,528 


0, 


,720 


2142.9 


3 


.190 


0, 


,371 


0, 


,712 


3 


,101 





,523 


0, 


,712 


2285.7 


2, 


.910 


0, 


,389 


0, 


,707 


2, 


,857 


0, 


,520 


0, 


,707 


2428.6 


2, 


.472 


0, 


,437 


0, 


,702 


2, 


,642 


0, 


,516 


0, 


,702 


2571.4 


2, 


.194 


0, 


,470 


0, 


,697 


2, 


,452 


0, 


,511 


0, 


,697 


2714.3 


2, 


.022 


0, 


,486 


0, 


,691 


2 


,282 


0, 


,505 


0, 


,691 


2857.1 


1, 


.905 


0, 


,499 


0, 


,685 


2, 


,133 


0, 


,513 


0, 


,685 


3000.0 


1, 


.818 


0, 


,506 


0, 


,678 


2, 


,031 


0, 


,515 


0, 


,678 


3776.8 


1, 


.527 


0, 


,498 


0, 


,646 


1, 


,610 


0, 


,498 


0, 


,646 


4754.7 


1, 


.199 


0, 


,502 


0, 


,624 


1, 


,221 


0, 


,494 


0, 


,624 


5985.8 


0, 


.909 


0, 


,491 


0, 


,597 


0, 


,880 


0, 


,489 


0, 


,597 


7535.8 


0, 


.667 


0, 


,481 


0, 


,563 


0, 


,630 


0, 


,484 


0, 


,563 


9487.0 


0, 


.440 


0, 


,500 


0, 


,545 


0, 


,430 


0, 


,493 


0, 


,545 


11943.5 


0, 


.304 


0, 


,473 


0, 


,533 


0, 


,272 


0, 


,475 


0, 


,533 


15036.0 





,210 


0, 


,457 


0, 


,511 





,166 


0, 


,465 


0, 


,511 


18929.2 


0, 


,145 


0, 


,448 


0, 


,480 


0, 


,111 


0, 


,439 


0, 


,480 


23830.6 


0, 


,100 


0, 


,424 


0, 


,445 


0, 


,075 


0, 


,417 


0, 


,445 


30001.0 


0, 


,069 


0, 


,400 


0, 


,420 


0, 


,033 


0, 


,400 


0, 


,420 
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scattered flux fraction, tlie direct attenuated stellar flux 
fraction, and the total escaping flux fraction, containing 
both attenuated stellar as well as scattered radiation, for 
the homogeneous case. Columns 7 to 10 repeat this in- 
formation for the equivalent clumpy case. All fluxes are 
referred to unity at each wavelength, with unit flux cor- 
responding to no attenuation, and in cases with clumpy 
dust distributions, they are averages over all directions. 

For a given geometry, the value of the V-band optical 
depth may be regarded as a measure of the dust mass 
present in that geometry. The attenuation optical depth 
(column 3 and 7) is defined as 



Tan = -ln(f(esc)) 



(1) 



and thus measures the absorption occurring in the sys- 
tem. This is a rather useful measurement since in galac- 
tic environments viewed from afar both the attenuated 
stellar flux as well as the scattered light are combined in 
the escaping total flux observed. Note that this defini- 
tion of the attenuation optical depth differ s from that of 



the effective optical depth (rcff) defined in WG96, which 
refers to the attenuated, non-scattered stellar flux only. 
The Tatt values are always smaller than the radial opti- 
cal depth listed in column 2, because the dust albedo is 
always non-zero and the stars are distributed throughout 
extended volumes. Also, the attenuation optical depth for 
the clumpy cases ne ver exc eeds that of the corresponding 
homogeneous cases ( WG96| ). It needs to be reemphasized 
that the clumpy case treated here is for one specific set of 
dumpiness parameters, namely a filling factor ff = 0.15, an 
interclump medium to clump density ratio k2/ki — 0.01, 
and a system size divided into 30"^ cubical bins. The varia- 
tions of transfer characteristics expected when any of these 
paramete rs are c hanged within reasonable ranges were ex- 
plored in WG9(: , and they are briefly discussed in Section 
3.5. 

3.2. Reddening Effects 

Reddening of a spectral energy distribution compared 
to that expected from spectroscopic indicators is usually 
taken as the first sign of the presence of dust in a galac- 
tic system. Color-color diagrams derived from multi-color 
photometry of galaxies frequently exhibit reddening ar- 
rows which purport to show the direction in which data 
points move in the color-color plane in response to the 
presence of a given column density of dust. Implicit in 
such displays is the assumption that the redde ning is pro- 
portional to the dust column. It was shown by WTC that 



the relationships between reddenings and dust columns is 
usually quite non-linear and highly dependent upon details 
of the geometry and the scattering properties of the dust. 
Here we show the additional effects due to the presence of 
a clumpy structure in the attenuating medium. 

In Fig. ^ we show the traditional reddening parameter 
E{B ~V) as a function of the homogeneous radial optical 
depth TV ; in Fig. ^ the congested region near the origin of 
Fig. ^ is enlarged. Several interesting points can be made. 
The most dramatic effect arising from the transition from 
a homogeneous structure to a two-phase clumpy structure 
occurs in the SHELL geometry. The homogeneous SHELL 
is similar to a spherical homogeneous screen, except that 
flux is scattered back into the beam in the SHELL geome- 
try. As expected for a screen, the predicted reddening for 



the homogeneous SHELL geometry increases linearly with 
increasing dust column density. The scaling factor is less 
than that for a screen due to the blue color of scattered 
radiation included in the measurement. Consequently, the 
ratios E{B — V)/tv for the homogeneous SHELL models 
are 0.20 for the MW case, 0.24 for the SMC case. This 
should be compared to the standard ratio of 0.35 for av- 
erage Galactic extinction. The transition to a two-phase 
clumpy dust structure in the shell leads to a reddening 
function which quickly becomes non-linear and saturates 
near E{B-V) - 0.25 for MW dust, near E{B-V) - 0.30 
for SMC dust. This initial saturation is caused by the fact 
that individual clumps become optically thick while the 
low-density interclump medium initially contributes little 
to the reddening. At Ty > 10, the clumpy SHELL geome- 
try exhibits a gradual increase in the reddening again, now 
solely due to the growing attenuation by the interclump 
medium. Note the interesting parallel to the behavior of 
a curve-of-growth in astronomical spectroscopy, where the 
increase of absorption in the core of a line is followed by 
saturation and subsequent further growth due to absorp- 
tion in the damping wings of the line. 

The DUSTY geometry exhibits a permanent saturation 
at larger optical depths both in the clumpy and in the 
homogeneous states, as was shown already by WTC] . The 
transition to a clumpy structure makes the saturation oc- 
cur at roughly half the amount of reddening in E{B — V) 
than was found in the homogeneous case. The saturation 
in the homogeneous case is a result of the fact that radia- 
tion emerges primarily from the outermost regions within 
optical depths one to two, regardless of how much dust is 
present further inside the system. This is fundamentally 
unchanged by adopting a clumpy structure, except that 
light now escapes preferentially through low-optical depth 
lines-of-sight between clumps, resulting in an overall lower 
amount of reddening. 

The level of E{B — V) at which realistic dusty struc- 
tures saturate is highly significant. It has been noted 
with great emphasis by Dickinson (1997), Sawicki and Yee 
(1998), and Steidel et al. (1999), among others, that the 
E{B — V) color excesses of Lyman-break galaxies crowd 
around values of E{B — V)^ 0.2, with hardly any exhibit- 
ing color excesses greater than 0.35. This result is fully 
consistent with our model results. All models except the 
homogeneous SHELL models saturate at reddening val- 
ues E{B — V) less than 0.3, a value which is breached 
by the clumpy SHELL models only for ry > 15. This 
suggests that the highly idealized homogeneous SHELL 
model, which is basically the historical screen model with 
scattering, does not apply to t hese L yman-break galax- 
ies, as was argued forcefully by WTC already. However, 
the result also means that the actual dust column densi- 
ties and the resulting attenuations are quite unconstrained 
by these color excess measurements. The fact that these 
galaxies were selected on the basis of their rest-frame-UV 
emission characteristics does argue, however, for moderate 
total attenuations. 

The saturation of color excesses occurs for higher col- 
umn densities of dust when wavelengths longer than the 
V-band are involved. In Fig. ||, this is illustrated in the 
case of a E{V~K) versus E{B — V) color-color plot. Here, 
models with increasing optical depths are plotted as series 
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of connected symbols, starting with Ty — at the ori- 
gin. The separation of symbols from one another along a 
given reddening line is a measure of how each color excess 
scales with optical depth. Not only is the reddening line 
not a simple arrow (with the exception of the homoge- 
neous SHELL model) , the detailed shape of the reddening 
line depends on the dust / star geometry and the structure 
of the dusty medium. Of particular interest again is the 
clumpy SHELL geometry, which shows the initial satu- 
ration in E{B — V) when clumps become optically thick 
in the V-band, followed by a saturation in E{V — K) at 
much higher radial dust column density when the same 
clumps reach the optically thick state in the K-band. Fi- 
nally, after clumps are sufhciently thick to block light in 
all bands, the interclump medium begins to dominate the 
reddening process. Where and when this occurs depends 
on the ratio of densities in the clump and the interclump 
medium. In a real physical system where the differences 
between different phases are rather more gradual than the 
fixed ratio assumed in our models, the transition shown in 
Fig. 5 between the reddening caused by the clumps and 
that caused by the interclump medium will also occur in 
a less abrupt manner. 



3.0 
2.5 
2.0 



i 1.5 



1 .0 



0.5 




0.0 

0.00 0.10 



0.20 0.30 
E(B-V) 
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Fig. 5.— E{V - K) is plotted versus E{B - V). The location of 
the symbols along each line indicated the change in colors as the op- 
tical depth is increased. For the SHELL geometry, the symbols are 
plotted for every model point. For the DUSTY and CLOUDY, 
the symbols are plotted for every fourth model Ty point. 

3.3. Attenuation Functions 

In Fig. o, we present the predicted attenuation opti- 
cal depths (see equation 1) as a function of inverse wave- 
length for MW dust for three cases, the optically thin 
case for Ty — 0.5 (top), the moderately optically thick 
case Ty — 1.5 (middle), and the very optically thick case 
Ty — 4.5 (bottom). Also shown in each panel is the av- 
erage galactic extinction for the same dust column densi- 
ties. While the left-hand panels show the absolute values 
of the respective optical depths, the corresponding right- 
hand panels present the same data normalized at V, in 
order to facilitate the comparisons of the changes in the 
wavelength dependence of each function. The DUSTY, 
SHELL, and CLOUDY models are contrasted in both their 
homogeneous and their clumpy forms. Fig. ^ presents the 



analogous set of data for SMC dust. Insert boxes are in- 
cluded to show the behavior of those curves which go out 
of bounds vertically on the main graphs. 

Several significant conclusions can be drawn from the 
illustrated cases. In all instances, attenuation optical 
depths for the three geometries are less than the inter- 
stellar extinction produced by the same column density of 
dust. This difference is due to the return of scattered 
radiation in all cases, exclusively so in the case of the 
SHELL geometry, and also due to lower attenuation in 
cases of stars only lightly embedded in the dust distribu- 
tion, which is happening in the CLOUDY and DUSTY ge- 
ometries. The descending order from SHELL, to DUSTY, 
to CLOUDY geometry can thus be understood. Simi- 
larly, clumpy structures containing the same amounts of 
dust as corresponding homogeneous structures always pro- 
vide less attenuation than the homogeneous ones at all 
wavelengths. More importantly, and this difference is par- 
ticularly prominent in the SHELL geometry, the clumpy 
structures provide a substantially grayer wavelength de- 
pendence for the attenuation than would have been ex- 
pected from adopting the original extinction curve. The 
trend toward gray attenuation increases strongly with in- 
creasing column density, because an increasing fraction of 
the total attenuation is caused by clumps which are opti- 
cally thick shortward of the V-band. For this same reason, 
the relative strength of the 2175 A UV extinction feature, 
prominent in the MW dust, is greatly reduced in clumpy 
structures. This is in contrast to the case of the homoge- 
neous SHELL model, which exhibits a 2175 A feature even 
stronger than that in the interstellar extinction (Fig. ^ 
RHS). This is due to the fact that the albedo of MW dust 
displays a minimum coinciding with this feature, leading 
to return less scattered light here than at wavelengths im- 
mediately longward or shortward of the feature. Similarly, 
the declining dust albedo shortward of 1500 A causes the 
predicted attenuation curve for the homogeneous SHELL 
model curves to be steeper in the far UV than the inter- 
stellar extinction curve. 

It is important to recognize that the slope of the attenu- 
ation functions in the UV for the more realistic embedded 
and clumpy geometries is dependent on the total dust col- 
umn density. The use of a single attenuation function, 
for example th e widely used "Calzetti Attenuation Law" 
(Calzetti 1997), for the dereddening and flux correction of 
individual starburst galaxies is therefore not justified in 
the absence of information about the dust column density. 
A comparison of the "Calzetti Attenuation Law" with pre- 
dictions from our models will be discussed in more detail 
in Section 4. 

3.4. Scattered Light Component 

The integrated light of dusty galaxies consists of a cer- 
tain fraction of scattered light, with a spectrum similar to 
that of the integrated star light. This fraction depends on 
wavelength, geometry, total dust column density, and dust 
scattering properties, as well as viewing direction in the 
case of disk galaxies. In the latter, the relative fraction of 
scattered light is largest in directions perpendicular to the 
disk, this being the direction in which the escape probabil- 
ity is highest for both stellar and particularly for scattered 
radiation. Given the spherical geometry of our models, we 
are able to predict the direction-averaged scattered light 
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Fig. 6. — The attenuation curves for the Milky Way model runs are plotted unnormalized (left) and normalized (right). 
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Fig. 7. — The attenuation curves for the Small Magellanic Cloud model runs are plotted unnormalized (left) and normalized (right). 
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Fig. 8. — The ratio of scattered to stellar flux (Fscat/^*) are plotted for the Milky Way (left) and SMC (right) model runs. 
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components. These are shown in Fig. ^ as fractions of the 
attenuated stellar flux for three characteristic visual radial 
optical depths (0.5, 1.5, 4.5), and for three geometries with 
both MW and SMC dust, both clumpy and homogeneous. 
As expected from observations of r eflection nebulae 



with embedded stars ( |Witt et al. 19924 1993), the scat- 
tered light fraction for homogeneous SHELL models easily 
exceeds unity and rises strongly with increasing optical 
depth, more as a result of the increasing attenuation of 
the starlight, less as a result of increasing scattered light 
flux. We view the embedded DUSTY, the clumpy SHELL 
models, and the CLOUDY models to be representative of 
actual geometries of galaxies. For these, the scattered light 
fraction appears to be fairly constant across the UV/visible 
spectrum and is likely between a quarter and one half of 
the total integrated hght. In the case of the MW dust, the 
depression near 4.6 /im~^ is due to the albedo minimum as- 
sociated with the 2175 A extinction bump, and this feature 
is actually seen in th e spectnmr of the diffuse (i.e. scat- 
tered) galactic light ( |Witt fc LiUie 1973| ). The transition 
from homogeneous to clumpy structures is associated with 
a substantial reduction of the fractional scattered light 
flux. This is due to two causes: Direct starlight escapes 
more easily from the system, mainly through gaps in the 
clump distribution, while the scattered light flux is dimin- 
ished, because clumps are mainly illuminated externally 
with the scattered light being directed preferentially into 
the clump interior where final absorption is more likely. 
In clumpy media, clumps in effect become the scattering 
centers instead of individual dust grains. As shown in 



WG96 , the effective albedo of clumps is always lower than 



the input dust albedo. Overall, the fractional amounts of 
scattered light for MW dust and for SMC dust are about 
the same for identical dust column densities. 

3.5. Variation in dumpiness Parameters 

We have emphasized that the models involving clumpy 
dust distributions presented in the paper are based upon 
one particular set of dumpiness parameters (f f — 0. 15, 
k2/ki = 0.01, N = 30). This set was found by |WG96| to 
be closest to that representing the structured ISM in the 
Galaxy. Detailed results on how the radiative transfer is 
affected when th e dumpiness parameters are varied can be 
found in WG96 as well. Here we summarize the principal 
effects. 

The density ratio k2/ki and the filling factor ff together 
control the contributions to the total opacity made by the 
interclump medium and the clumps. Our case oi ff = 0.15 
and k2/ki — 0.01 is intermediate to the case of ^2/^1 near 
unity, where the interclump medium controls the system 
opacity, with results very similar to the homogeneous case 
(^2/^1 = 1-0) and the case of fc2/^i < 10^"^, where the 
opacity is controlled by the blocking effect of the clumps 
and is then proportional to the filling factor. As the fill- 
ing factor increases, the appearance of the system will ap- 
proach that of the homogeneous case. 

The grayness of the attenuation (Section 3.3) increases 
as k2/ki decreases for the same amount of dust. Also, the 
ratio of scattered light to direct, partially attenuated stel- 
lar light (Fig. 8) decreases as the ratio k2/ki decreases, 
both as a result of the growing probability of stellar light 
escaping without any interaction and as a consequence of 
the lower scattering ability of a medium dominated by 



increasingly optically thick clumps. Finally, for a given 
amount of dust, the overall opacity of the medium will 
increase as the number of clumps (N) increases. 

4. DISCUSSION 

In this section we discuss which of our models with re- 
spect to geometry, optical depth, and dust type might best 
be used for the analysis of observations of different kinds 
of galaxies. In particular, we concentrate on dust attenua- 
tion studies of nearby starburst galaxies and Lyman break 
galaxies. 

4.1. The "Calzetti Attenuation Law" 

Starburst galaxies exhibit a wide range of dustiness as 
evidenced by large variations in the UV spectral index (3, 
ratios of recombination line strengths, and ratios of far-IR 
to UV luminosities. The UV spectral index (3 was defined 
by Calzetti et al. (1994) through the relation F{\) oc A'' in 
the range 1250 < A < 2600 A. By comparing the spectral 
energy distributions (SED) of more highly reddened star- 
bursts with those of barely reddened starbursts, Calzetti 
et al. (1994) derived the observed wavelength dependence 
of the dust attenuation in these systems. The attenua- 
tion function takes into account the return of scattered 
light into the integrated light of a galaxy and also reflects 
radiative transfer effects arising in the three-dimensional 
distributions of stars, gas, and dust in these galaxies. This 
observed attenuation function, further refined in Calzetti 
(1997), is referred to in the literature as the "Calzetti At- 
tenuation Law" or "Calzetti Law". This attenuation due 
to internal dust must be clearly distinguished from the 
wavelength dependence of interstellar extinction, which 
is the attenuation due to absorption and scattering by a 
screen of dust located between a distant source and an 
observer. 

The principal characteristic of the "Calzetti Law" is the 
absence of a 2175 A feature. Gordon et al. (1997) showed 
that the "Calzetti Law" can be reproduced, if the radia- 
tive transfer in starburst galaxies is occurring in a medium 
with intrinsic dust properties closely resembling that of 
SMC dust. Also, the large range of observed reddenings, 
as measured by variations in the spectral index (3, appears 
to require efficient attenuation geometries such as provided 
by our SHELL or DUSTY geometries. 

In Fig. ^ we compare attenuation values derived from 
starburst galaxies by Calzetti (1997) with our model at- 
tenuation functions. In order to facilitate this com- 
parison, models and observations are renormalized to 
E{\ — V)/E{B — V). In order to transform the values 
in Table 5 of Calzetti (1997), we have used the result 
of Calzetti (1997) that the stellar continuum suffers 40% 
the reddening that the gas suffers, i.e. E{B — T^)stciiar = 
QAE{B — V)gas- The data point capable of discriminat- 
ing between different models most clearly is the one near 
6 /^m~^, derived from the far-UV slope in the reddened 
SEDs of starbursts. No model with MW dust matches this 
data point. In addition, the lack of a 2175 A feature in the 
observations compared to the prediction of a 2175 A fea- 
ture in the models with MW dust argues strongly against 
the presence of MW dust in starburst galaxies. This will 
be discussed further in Section 



4.4 



Among SMC dust models, the same point does discrim- 
inate between dust column densities and clearly favors the 
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Fig. 9.— The _E{A - V) /E(B - V) curves for MW and SMC arc plotted along with the Calzetti (1997) data (filled squares) and the Calzetti 
et al. (1994) data (solid line) both derived from starburst galaxy observations. 
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clumpy SHELL model for Ty — 1-5. As is apparent from 
Fig. I?! such a model predicts an attenuation optical depth 
of about Tatt = 2.0 at 6 ^m~^, and, according to Fig. |[ 
a color excess E{B — V) ~ 0.22. An attenuation optical 
depth of 2.0 implies that the flux at 6 fim~^ has been re- 
duced by a factor 7.4 on average in this UV-selected sample 
of starburst galaxies . We need to re-emphasize that there 
is no universal attenuation function and that the atten- 
uation factor of 7.4 corresponds simply to the average of 
the sample of UV/optically-selected starbursts studied by 
Calzetti et al. (1994). Our average attenuation factor of 
7.4 compares well with the luminosity-weighted mean ab- 
sorption factor of 5.4 derived by Meurer et al. (1999) for 
the same sample. We discuss the UV flux reduction fac- 
tor and the uncertainties associated with its determination 



further in Section 4.2 
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Fig. 10.— The ratio of FIR to UV-NIR flux (Ffir/Fuv-nir) is 
plotted versus E{B — V). Each model curve is shown twice, once as- 
suming a 1 million year old constant star formation SED and once 
assuming a 1 billion year old constant star formation SED. The 
symbols are plotted every third model ry point. 

4.2. The Far-IR to Optical Flux Ratio 

It came as an initial surprise during the IRAS mission 
that far-IR ultraluminous and by inference, very dust-rich, 
galaxies in most instances exhibited optical colors which 
were barely distinguishabl e from those of normal galaxies 
of similar population type. WTC showed that this was due 



to the fact that, given embedded dust geometries, the op- 
tical colors were determined by the least attenuated stars, 
while the far-IR flux was due to the absorption of light 
from stars which contribute little or nothing to the optical 
flux. The present calculations were extended to dust col- 
umn densit ies cor responding to ry = 50 instead of 15 in 
the case of WTC, In Fig. |lO[ we show the resulting ratios 
of the far-IR flux, integrated over all frequencies at which 
dust emission occurs, to the integrated UV/optical/near- 
IR stellar flux. We reproduce the most extreme observed 
ratios near 100 with the clumpy SHELL model of ry = 50 
with a maximum reddening near 0.8 in E{B — V) and only 
slightly lower ratios with the homogeneous DUSTY mod- 
els with a reddening less than 0.3 in E{B — V). Only the 
homogeneous SHELL, as expected for a simple screen, pro- 
duces very large amounts of reddening of E{B — V)^ 4.0 
and FiB./Fuv-NiR ~ 1000 for ry = 50. The fact that 



most IR luminous galaxies appear to be relatively blue 
seems to indicate that clumpy dust structures prevail in 
most of them. This is also confirmed by the correlation 
between the ratio Fir/Ffi60BW and the spectral index Ad 
for different geometries and dust types, shown in Fig. n3. 

4.3. The UV Flux Attenuation Factor 

A major aim of studies of high-z starbursts (or Lyman- 
break galaxies) is to determine the epoch of peak star 
formation in the history of the Universe, if such a peak 
indeed exists. The outcome of such an investigation de- 
pends critically upon the corrections applied for dust at- 
tenuation in the rest-frame UV continuum spectrum of 
the observed galaxies, given that these systems contain 



substantial amoun t s of internal d u st ( Dickinson 1997 



Meurer et al. 1997 ; Heckman 199g ; Meurer, Hcckman, & 
Calzetti 1999| ; [Stcidcl ct al. 19991 ). In addition, the upcom- 
ing NASA GALEX mission is to produce about 100,000 
UV galaxy spectra in order to probe the global history of 
star formation, and the observed UV spectral slope is be 
be related uniquely to the attenuation factor by which the 
UV flux has been reduced by dust internal to the galaxies. 

Our Figures ^ and ^ already indicated some of the prob- 
lems to be expected in such an approach. The attenuation 
functions are geometry dependent and they become in- 
creasingly gray as the dust column density increases. This 
means that a heavily attenuated starburst may exhibit a 
steep UV slope in its SED, resembling that of a barely 
reddened galaxy, but show a greatly reduced UV flux, cor- 
responding to a system with little star forming activity. 
The chances for underestimating the star formation rate 
in early epochs of the Universe are therefore severe, if dust 
is abundant generally in the early galaxies. 

Fortunately, our models provide an excellent set of tools 
with which to evaluate different approaches to the problem 
of determining the far-UV attenuation factor. In Fig. jlja 
we illustrate the computed attenuation factor at 1600 A~as 
a function of E{B — V) found for different model geome- 
tries and dust types. As anticipated above, the observation 
that most observed values of E{B — V) for Lyman-break 
galaxies fall into the narrow range of 0.2 to 0.3 barely con- 
strains the attenuation factor to the range 4 to 200. This is 
not very useful as a constraint. In Fig. pTp we investigate 
the use of the spectral index /3 as a tool to evaluate the 
attenuation factor. Our radiative transfer models do not 
assume a specific energy distribution for the source spec- 
trum. Hence, we predict only the change in the spectral 
index /3 which occurs as a result of the wavelength depen- 
dent attenuation. The GALEX program advocates this 
particular approach by suggesting that the change in the 
UV slope of the galaxies' SED can be converted into infor- 
mation about the applicable attenuation factor. Meurer 
et al. (1999) have taken the additional step of calibrat- 
ing the correlation between the spectral index /3 and the 
UV-attenuation, derived for nearby UV-selected starburst 
galaxies. In observed starburst spectra, this change in 
spectral index ranges from zero to 3 (Meurer et al. 1997). 
Our Fig. pT| b reveals that correlations between the spectral 
index change and the UV attenuation factor are extremely 
model dependent and even more dependent on the type of 
dust present in a given system. The calibration proposed 
by Meurer et al. (1999) is shown in Fig. pT^ j as a solid black 
line. For this to be valid universally, one must assume that 
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Fig. 11.— The attenuation factor in the HST/WFPC2 F160BW filter is plotted versus E{B - V) (a) and A/3 (b). The symbols are plotted 
every third model Ty point. 



SMC-type dust is present in all starburst galaxies and that 
SHELL-type geometries with various degrees of dumpi- 
ness are the prevailing star/dust geometries. Even when 
restricting considerations to SMC-type dust, the range of 
observed values of A/3 is consistent with attenuation fac- 
tors in the range 0-150. As is evident from Fig. [ii}d, MW- 
type dust is not expected to produce positive values of 
A/3. This is contrary to observations in starbursts, where 
such positive changes in the UV spectral index are seen 
(Meurer et al. 1999). However, MW-type dust will reveal 
its presence by the appearance of a 2175 A feature. This 
has bee n seen in the heavily obscured starburst nucleus 
of M33 ( [Gordon ct al 1999a[ ), which indicates that SMC- 
type dust is not a universal characteristic of all starburst 
environments. However, the preference for SMC-type dust 
in starburst galaxies is affirmed, when we compare values 
of -Ffir/^fi60BW flux ratios and corresponding values of 
A/3 from observations (Meurer et al. 1997) with our mod- 
els in Fig. O. For this calculation we assumed that the 
unattenuated source SED was given by a 1 or a 1000 Myr 
constant star formation SED of Fioc & Rocca-Volmerange 
(1999). The observed data points clearly follow the trend 
set by the model for clumpy SHELL configurations with 
SMC dust. This agreement would become even more con- 
vincing, were we to present models with different ratios of 
clump to interclump density, different clump size, and dif- 
ferent filling factors, which is, however, beyond the scope 
of this present discussion. The demonstrated model de- 
pendence of the ratio of FIR to F160BW flux suggests 
that the observed scatter in the data points is most likely 
real, as was concluded by Meurer et al. (1999) also. The 
spread in the FIR to F160BW ratio for a given value of 
A/3 corresponds to a factor of 10 uncertainty in the cor- 
responding UV attenuation factor. We conclude, therfore, 
that the UV spectral slope is not a reliable indicator of the 
dust attenuation in a given galaxy. It is usable only, if a 
large sample of starburst galaxies with relatively uniform 
characteristics is being considered at once. 

In Fig. p^a, we plot the ratio of dust-emitted IR flux 
to the integrated UV-NIR flux against the attenuation 
factor at 1600 A. The model calculations assume two star- 
bursts of different ages, one a 1 million year old constant 



star formation SED, the other a 1 billion year old con- 
stant star formation SED. We note that even this ratio 
does not yield a unique relation to the UV attenuation 
factor. In particular, dumpiness in a given dust distribu- 
tion will reduce the attenuation at 1600 A for a given value 
of the ratio of dust-emitted IR flux to the integrated UV- 
NIR flux, because dumpiness is associated with a much 
grayer attenuation function. This is particularly evident 
in SHELL-type geometries. Variations in the age of the 
starburst and the associated changes in the UV-NIR SED 
affect the ratio as well, leading to a significant model de- 
pendence of the predicted attenuation factor at 1600 A. 
By contrast, if we plot the ratio of integrated FIR flux to 
the measured flux at 1600 A against the attenuation factor 
at 1600 A, as shown in Fig. [l2|b, the model and age de- 
pendence practically vanish. We conclude, therefore, that 
the UV attenuation factor can be evaluated reliably for 
individual galaxies only, if a complete spectral survey of 
the FIR emission of the galaxy is carried out and a flux 
measurement in the far-UV is obtained. 
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Fig. 12. — The far-IR to UV-NIR (a) and far-IR, to F160W (b) flux ratios are plotted versus the attenuation factor at 1600 A. Each model 
curve is shown twice, once assuming a 1 million year old constant star formation SED and once assuming a 1 billion year old constant star 
formation SED. The symbols are plotted every third model Ty point. 



Fig. 13.— The ratio of FIR to F160BW flux (FpiR/FpigoBw) is 
plotted versus A/3. Each model curve is shown twice, once assuming 
a 1 million year old constant star formation SED and once assum- 
ing a 1 billion year old constant star formation SED. The symbols 
are plotted every third model Ty point. The solid black points are 
from Meurer et al. (1999) and i^pigoBW ~ ^1600 = ^cfieoo and 
Ac = 1600 A. The A/3 values for the galaxies were computed using 
A/3 = /3 + 2.5. 

4.4. Lack of 2175 A Extinction Bump in Starbursts 

Gordon et al. (1997) showed that the absence of the 
2175 A extinction bump in the attenuation function of 
starbursts cannot be explained by radiative transfer effects 
but requires the actual absence of the absorber giving rise 
to this feature. This is despite the fact that some weaken- 
ing of the 2175 A feature was seen in the attenuation curves 
(Fig. H) for embedded geometries involving MW dust and 
further weakening occurred when dumpiness was included. 
As was found by Gordon et al. (1997), the 2175 A bump 
can be suppressed using radiative transfer effects, but at 
the expense of a gray extinction throughout the UV. Thus, 
the existence of a substantial slope in the Calzetti (1997) 
attenuation curve (Fig. ||) implies that the dust in star- 
burst galaxies must lack the absorbers which produce the 
2175 A bump. The dust in sta rbursts most closely resem - 
bles that found in the SMC ( [Gordon fc Clayton 199^ . 
Recent modeling of this dust by Zubko (1999) shows that 
such dust lacks the component normally introduced to pro- 
duce the 2175 A feature (graphite), and that the overall 
size distribution of SMC grains is shifted toward smaller 
sizes compared to MW dust. 

We find the color index (F170M^ - F218W) to be par- 
ticularly sensitive for distinguishing between the presence 
of SMC dust and MW dust, one without the 2175 A fea- 
ture and the other with. In Fig. we show a color-color 
diagram of starburst data, which are compared to the col- 
ors of unreddened starburst models (open symbols) and to 
reddening lines produced by our models. These reddening 
lines are attached to different unreddened models to show 
them distinctly; they may be shifted horizontally to fit the 
data. Several facts emerge clearly: All starburst data ap- 
pear on or below the line of unreddened burst models, i.e. 
they are reddened in the {F170W - F218M^) index. Only 



the reddening lines for SMC dust fall below the unred- 
dened bursts, while all MW models fall above the unred- 
dened bursts, where no data points are found. Both the 
amount of reddening in the (F170I^ — F218W) index and 
the distribution of data points in the color plane can best 
be fitted with a clumpy SHELL model with SMC dust. 
We take this (as did Gordon et al. 1997) to be one of the 
strongest indications that nearby starbursts contain SMC 
type dust and that the prevailing geometry is that of a 
clumpy SHELL model. This sensitive test might be ap- 
plied to the rest-frame UV SEDs of high-redshift galaxies 
in general in order to test for the presence of the type dust 
(MW or SMC). 
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Fig. 14.— The (F170W-F218W) and (B-V) colors for starburst 
observations and constant star formation SES models are plotted 
versus each other along with the reddening trajectories predicted 
from our radiative transfer model. Each reddening trajectory could 
be attached to any of the SES model points. For this figure, we 
have arbitrarily attached placed the reddening trajectories to define 
the envelope of the starburst observations. The starburst observa- 
tions all lie above [(F170W-F218W) > -0.6] the SES models, which 
is only possible if the dust in starburst galaxies is similar to that 
found in the SMC. 

5. SUMMARY AND CONCLUSIONS 
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We have presented new multiple-scattering radiative 
transfer calculations for three representative types of 
galactic environments and two dust types, as found within 
the Local Group in the form of the average MW and SMC 
type dust. Both homogeneous and two-phase clumpy dust 
distributions have been explored over a large range of 
dust column densities, measured by ry = — 50. The 
full numerical details of the models are made available 
to the community via the authors, the authors WWW 
site, and eventually various national/international data 
centers. The environments examined include DUSTY, a 
fully mixed spherical distribution of dust and embedded 
stars; SHELL, a dust-free star filled region, surrounded 
by a dust shell; and CLOUDY, a mixed dust-star region 
surrounded by stars in a dust-free region on the outside. 
The models cover the wavelength range 1000 A to 30000 
A with exceptionally close coverage of the rest-frame UV. 
This aspect makes these models particularly applicable to 
UV-selected starburst galaxies and Lyman-break galaxies. 

We examined the models with respect to predicted red- 
dening effects, attenuation functions, and the contribution 
of scattered light to the total flux as a function of wave- 
length. All model types except the homogeneous SHELL 
model exhibit reddening saturation effects. In particular, 
models that are likely representatives for starburst galax- 
ies show values of E{B — V) not in excess of 0.3, despite a 
large range of dust column density. Reddening "vectors" 
in color-color diagrams are complex non-linear functions 
depending on the prevailing dust type and structure, as 
well as on the geometry. Attenuation functions, which 
measure the wavelength-dependent reduction in the to- 
tal stellar and scattered flux from a dusty galaxy, exhibit 
a trend toward increasing grayness with increasing dust 
column density. This trend is stronger in clumpy dust 
distributions than in homogeneous ones. No justification 



was found for the use of a universal attenuation function 
for the analysis of a large sample of galaxies. In particu- 
lar, we found that the widely-employed "Calzetti Law" is 
most closely reproduced by our clumpy SHELL model with 
SMC-type dust and a dust column density corresponding 
to Ty = 1.5. Such a case corresponds to an attenuation 
optical depth in the UV of Tatt — 2.0 or a UV attenuation 
correction factor of 7.4. The Calzetti law is valid for sta- 
tistical samples which have similar dust/star geometries 
as found in UV-selected starbursts and which have dust 
column densities within narrow limits. Under no circum- 
stances should it be used for the attenuation correction for 
single galaxies. 

When investigating various approaches to determining 
this attenuation factor in general, we found that the use of 
the color excess E{B — V) and the UV spectral index f3 lead 
to results that depend strongly on dust type, dust distribu- 
tion, and source-dust geometry. Only the measurement of 
the -FfirZ-Ffigobw flux ratio promises reasonable certainty 
for the determination of the UV attenuation correction 
factor in individual galaxies. This reaffirms the need for 
observations of starburst galaxies and Lyman-break galax- 
ies over the entire range of rest-wavelengths, particularly 
the FIR, in addition to the UV. We further showed that 
observations can clearly discriminate between the presence 
of MW-type dust and SMC-type dust through the use of a 
(i^l70M^-F218W^) vs. {B-V) color-color diagram. When 
applied to existing data on UV-selected nearby starburst 
galaxies, a consistent prevalence of SMC dust is found. 
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Geoffrey Clayton, and Karl Misselt are also gratefully ac- 
knowledged. 



REFERENCES 



Berlind, A. A., Quillen, A. C, Pogge, R. W., & Sellgren, K. 1997, 
AJ, 114, 107 

Bianchi, L., Clayton, G. C, Bohlin, R. C, Hutchings, J. B., & 

Massey, P. 1996, ApJ, 471, 203 
Boisse, P. 1990, A&A, 228, 483 

Bruzual, A. C, Magris, G., & Calvet, N. 1988, ApJ, 333, 673 
Byun, Y. I. Freeman, K. C, & Kylafis, N. D. 1994, ApJ432, 114 
Calzetti, D. 1997, AJ, 113,162 

Calzetti, D., Bohlin, R. C, Gordon, K. D., Witt, A. N., Bianchi, L. 

1995, ApJ, 446, L97 
Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 

582 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
[CCM] 

Corradi, R. L. M., Beckman, J. E., & Simonneau, E. 1996, MNRAS, 
282, 1005 

di Bartolomeo, A., Barbaro, G., & Perinotto, M. 1996, in New Light 
on Galaxy Evolution, lAU Symp. 171, ed. R. Bender and R. L. 
Davies, Kluwer Acad. Publ., 363 
Dickinson, M. 1997, in The Hubble Deep Field, STScI Symp., eds. 

M. Livio, S. M. Fall, & P. Madau, (in press) 
Ferrara, A., Bianchi, S., Cimfitti, A., &z Giovanardi, C. 1999, ApJS, 



in press ( |astro-ph/9903078[ ) 
Fioc, M. & itocca- Volmerange, B., 1999, in preparation 
Fitzpatrick, E. L. 1999, PASP, 111,63 



Gordoi L, K. P., Calzetti, D., fc Witt, A. N. 1997,ApJ, 487, 625 
Gordon, K. D. k Clayton, G. C. 1998, ApJ, bUU, 81b 
Gordon, K. D., Witt, A. N., Carruthers, G. R., Christensen, S. A., 

& Dohne, B. C. 1994, ApJ, 432, 641 
Gordon, K. D., Hanson, M. M., Clayton, G. C, Rieke, G. H., & 

Misselt, K. A. 1999a, ApJ, 519, 165 
Gordon, K. D., Witt, A. N., & Clayton, G. C. 1999b, in preparation 
Heckman, T. M. 1998, in The Birth of Galaxies: Xth Recontres de 
Blois, ed. B. Guiderdoni, in press 



Henyey, L.G. & Greenstein, J.L. 1941, ApJ, 93, 70 
Hobson, M. P. & Scheuer, P. A. G. 1993, MNRAS, 264, 145 
Holmberg, E. 1958, Medd. Lund. Astron. Obs. Ser. 2, No. 6 
Holmberg, E. 1975, Stars and Stellar Systems IX, Univ. of Chicago 

Press 123 
Hurwitz, M. 1994, ApJ, 433, 149 

Kim, S.-H., Martin, P. G., & Hendry P. D. 1994, ApJ, 422, 164 
Kuchinski, L. E., Terndrup, D. M., Gordon, K. D., & Witt, A. N. 

1998, AJ, 115, 1438 
Larson, R. B. & Tinsley B. M. 1978, ApJ, 219, 46 
Lehtinen, K. & Mattila, K. 1996, A&A, 309, 570 
Lillie, C. F., & Witt A. N. 1976, ApJ, 208, 64 

Meurer, G. R., Heckman, T. M., Lehnert, M. D., Leitherer, C, & 

Lowenthal, J. 1997, AJ, 114, 54 
Mquror, G. R., Hock,man, T. M., & Calzetti, D. 1999, ApJ, in press 



( astro-ph/9903054 ) 
Misselt, K. A., Clayton, G. C, & Gordon, K. D. 1999, ApJ, 515, 128 
Morgan, D. H., Nandy K., & Thompson G. L 1976, MNRAS, 177, 



531 
Natta, A. i 
Sawicki, M 
Sawicki, M 

Smith, D. A., neiLei, CtC iicLyiies, ivi. r i xijfu^o, J^pj , 

Steidel, C. C. et al. 1999, ApJ, in press (^stro-ph/9811399| ) 
Takagi, T., Arimoto, N., & Vansevicius, V . 1999, ApJ, m press (astro- 



Panagia, N. 1984, ApJ, 287, 228 
& Yee, H. K. C. 1998, AJ, 115, 1329 
J., Lin, H., & Yee, H. K. C. 1997, AJ, 113, 1 
Herter, T., & Haynes, M. P. 199S, ApJ, 494, 150 



ph/9902219 ) 

Toller, G. N. 1981, PhD Thesis, State Univ. of New York at Stony 



Dwek, E. 1999, ApJ, submitted 



Brook 
Varosi, F.. 

Witt, A. N. 1977, ApJS, 35, 1 
Witt, A. N., Friedmann, B. C, Sasseen, T. P. 1997, ApJ, 481, 809 
Witt, A. N. & Gordon, K. D. 1996, ApJ, 463, 681 [WG961 
Witt, A. N. & Lillie, C. F. 1973, A&A, 25, 397 



18 



Clumpy Dust. II. 



Witt, A. N., Lindell, R. S., Block, D. L., & Evans, R. 1994, ApJ, 
427 227 

Witt, A. N., Oliveri, M. V., & Schild, R. E. 1990, AJ, 99, 888 
Witt, A. N., Petersohn, J. K., Bohlin, R. C, O'Connell, R. W., 

Roberts, M. S., Smith, A. M., & Stecher, T. R 1992a, ApJ, 395, 

L5 

Witt, A. N., Petersohn, J. K., Holbcrg, J. B., Murthy, J., Dring, A., 
& Henry, R. C. 1993, ApJ, 410, 714 



Witt, A. N., Thronson, H. A., & Capuano, J. M. 1992b, ApJ, 393, 
611 [WTC] 

Witt, A. N., Walker, G. A. H., Bohlin, R. C, & Stecher, T. P. 1982, 

ApJ, 261, 492 

Wolf, S., Fischer, O., & Pfau, W. 1998, A&A, 340, 103 
Zubko, V. G. 1999, ApJ, in press 



